Chemical equilibration and thermal dilepton production from the 
quark gluon plasma at finite baryon density 



OS 



in 



X 



D. Dutta, K. Kumar, A. K. Mohanty and R. K. Choudhury 

Nuclear Physics Division, Bhabha Atomic Research Centre, Bombay-400085 

(February 1, 2008) 



Abstract 

The chemical equilibration of a highly unsaturated quark-gluon plasma has 

been studied at finite baryon density. It is found that in the presence of small 

amount of baryon density, the chemical equilibration for gluon becomes slower 

and the temperature decreases less steeply as compared to the baryon free 

plasma. As a result, the space time integrated yield of dilepton is enhanced 

l/"") ' if the initial temperature of the plasma is held fixed. Even at a fixed initial 

Q^ ■ energy density, the suppression of the dilepton yields at higher baryo-chemical 

OS ■ 

potential is compensated, to a large extent, by the slow cooling of the plasma. 
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qj\ I. INTRODUCTION 



One of the important objectives of all the future collider experiments at RHIC and LHC 
is to detect the new state of matter called quark gluon plasma (QGP) which is expected 
to be produced during the ultra-relativistic heavy ion collisions. Such an exotic state, if 
formed, will cool during expansion till it reaches a critical temperature T c where the QGP 
will be transformed to the hadron phase via a first or second-order phase transition. The 
hadrons will then cool down further to the freeze-out temperature Tf where they cease to 
interact with each other and fly away to the detectors. During the process of the thermal 
expansion, photons and dileptons are produced directly from the plasma as well as from the 
hadron phases. These thermal photons and dileptons are considered to be ideal probes for 
the detection and study of subsequent evolution of the QGP [JTJ] as they leave the plasma 
without any interaction. 

In the standard scenario @ , the quark gluon plasma formed during the collision is expected 
to thermalize in a typical time scale of 1 fm/c. However, some recent works ||^[] suggest 
that due to high initial parton density (mostly gluonic) at RHIC and LHC energies, the 
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plasma may attain kinetic equilibrium in a very short time (r pa 0.3/m — 0.7 fm) but it 
may be far from chemical equilibrium. Since the initial parton plasma is mostly gluonic, 
many more quarks and anti-quarks are needed in order to achieve chemical equilibrium. 
Earlier studies have shown that a chemically non-equilibrated plasma cools faster as 
compared to equilibrated plasma, which follows Bjorken's scaling law (T 3 r=const ). In the 
case of non-equilibrated plasma, additional energy is consumed in producing more quarks 
and anti-quarks in approaching chemical equilibrium due to which the cooling is accelerated 
as compared to the Bjorken's scaling law. Since the plasma cools faster, it may not attain 
chemical equilibrium before the temperature reaches the critical value T c . Even if one 
includes transverse expansion, the large velocity gradient may drive the system further 
away from the chemical equilibrium as shown by ||. 

The effect of chemical equilibration on thermal photon and dilepton emission has been 
studied by several authors f7HT3j. In these studies, it was assumed that the nucleus- nucleus 



collision is fully transparent, leading to the formation of a baryon-free plasma. However, 
some of the recent works have suggested that even at RHIC energies the colliding nuclei 
may not be fully transparent and some amount of baryon stopping may occur particularly 
at higher rapidities . For Au + Au collisions at 200 A GeV, the baryo-chemical potential 



Hb may be of the order of T for y pa 2 |L4]], which is the region of interest of the PHENIX 
detector for di-muon measurements [15]. The presence of finite baryon density may affect 
the process of chemical equilibration and the rate of cooling. This will, in turn, affect the 
thermal dilepton yields, which are the potential probes for the detection of the expanding 
quark gluon plasma. 

The present work is aimed at studying these effects, by considering the dynamical evo- 
lution of a baryon rich plasma undergoing chemical equilibration. It is seen that the rate of 
cooling of the plasma slows down in the presence of /is compared to a baryon-free plasma. 
The calculations also show that the rate of dilepton production is suppressed at finite baryon 
density, but the time integrated thermal yields are rather unaffected due to slower cooling 
of the plasma in the presence of small amount of the baryon density. 

II. CHEMICAL EQUILIBRATION 

We assume a quark gluon plasma which has thermally equilibrated but is far off from 
chemical equilibrium. The distribution functions for quarks (anti-quarks) and gluons for 
such an unsaturated plasma can be described by Juttner distributions given by 
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f (_ \ _ _ \{g) e±X . f -^g / 1 \ 

M?; A g ( g -) + e^)/T A5( . )e ±x + eP /T > /a e p/T — X g {) 

where x = njT. The fugacity factor A, (i — q, q and g) gives the measure of the deviation 
of the distribution functions from the equilibrium values and \i (= /is/3) is the quark- 
chemical potential. The chemical equilibrium is said to be achieved when Aj — ► 1. It may 
be mentioned here that one can also define the quark and anti-quark distribution functions 
using a different definition for fugacities Xq and Xq given by 

X Q(Q) = e±X \{q) ( 2 ) 

In that case, one does not use the quark chemical potential /i explicitly. However, the 
definition (1) is quite convenient, since at equilibrium (A 9 =Ag=l), the chemical potential 
associated with Aj vanishes but the baryo chemical potential still exists. Further, we have 
taken X q = Xg at all values of r so that when \i — > 0, Xq = Xq = X q resulting in a baryon 
symmetric matter. 

A. Rate equations 

In general, chemical reactions among partons can be quite complicated because of the 
possibility of initial and final state gluon radiations. However, we restrict our considerations 
to the following dominant reaction mechanisms for the equilibration of parton flavours, 

99 ^ 999 99 ^ QQ 
The evolution of the parton densities are governed by the master equations 

3 p (n/) = (#2^3 - i2s- 2 ) " {Rg^g ~ Rg^g) (3) 

d,(n q un = d„(n q u») = (R g ^ q - R q ^ g ) (4) 



t-g-^q 



where i?2^3 and i?3^2 denote the rates for the process gg — > ggg and its reverse and R g 
and R q ^ g for the process gg — > qq and its reverse respectively. In case of a baryon rich 
plasma, the presence of quark-chemical potential //, affects the process gg ^ qq directly. 
We, therefore, study the \i dependence of this process more explicitly. The rate equations 
for this process are given by |16[ 



R °-« = 2J JW2eJ JW2eJ JW2eJ 72^m (27r) 5 (Pi+^-ft-^) X 
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E \M„ q \ 2 f g ( Pl )f g (p 2 )[l - f q (p 3 )][l - fM] (5) 

and 

1 f d 3 pi f d 3 p 2 f d 3 p 3 f d 3 p 4 , .44 



R ^ = 2 J T2^2eJ wtieJ wnKj wriE^ 6 x 

£ \M q ^ g ff q (p 3 )fg( P4 )[l + f g ( Pl )][l + f g ( P2 )} (6) 

In (5), the squared matrix element, summed over spin and color, J2 \M\ 2 is weighted by two 
gluon distribution functions f g for the initial states. The factor [1 — f q \ [1 — f q ] indicates 
Pauli blocking for the final states. In the reverse process (6), the rate is weighted by the 
distribution functions of quarks and anti-quarks for the initial states and the gluon final 
states each gain an enhancement factor [1 + f g ] due to Bose-Einstein statistics. The factor 
of 1/2 accounts for the identity of the two gluons. 
Using the identity 

[i - /,] [i - M = AMM^il e fe +P4 )/T (7) 

\A q 

and 

[i + fM][i + fM] = AM^M e ( P1+P2 )/T (8) 

\ 

it may be shown that the \x dependence in the rate equations (5) and (6) for the process 
99 ^ <?<? basically comes through the product of quark and anti-quark distribution functions 

fqfq- 

Fig. 1(a) and (b) show the variation of the product f q f q with momenta p and \ q respec- 
tively for different values of x. Fig. 1(a) shows the product f q f q as a function of p at a few 
typical values of x (= n/T) for X q — 0.1 and for temperatures T=0.57 GeV and T=0.2 GeV. 
Fig. 1(b) shows the f q f q as a function of X q at temperature T= 0.57 GeV and p=0. It is 
seen (Fig la) that as the momentum increases or the temperature decreases, the e P//T factor 
in the denominator of (1) dominates and the product f q f q becomes less sensitive to fi. Here, 
we have used the same values of p for both quark and anti-quark distribution functions, al- 
though their momenta p q and p q can be quite different. However, for any other combinations 
of p q and p q , the deviation will not be more than that observed for p=0. We have taken 
Pq = Pg = i n Fig. 1(b) to show the maximum amount of deviation of this product as a 
function of X q . It can be seen from Fig 1(b), that the above product is not very sensitive 
quark chemical potential if the plasma is highly unsaturated (X q < 0.1). This feature has 
an important relevance for the collisions at RHIC energy as the chemical equilibration is 



4 



not achieved by the time T drops to T c and the quark fugacity remains much below unity 
(X g ~0.1) ||. The nearly x independence of the above product is also evident from eq. (1). 
For small values of X g , the contribution from the factor X q e ±x in the denominator is not very 
significant if x is small. The x dependence of the quark and anti-quark distribution functions 
mainly arises due to the e ±x factor in the numerator. Since these exponential factors get 
cancelled in the product, f q f q will have weak x dependence at small baryon density if the 
plasma is highly unsaturated. So the integrals (5) and (6) will not depend explicitly on /i 
except through the thermal quark mass m 2 which is generally used as a cut-off parameter in 
the evaluation of the integral to avoid the divergence. The thermal quark mass appropriate 
for a non-equilibrium plasma is given by []T 7 7| , p~8 

,2 
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j o dpp[2 f g + f g + M (9) 



Since the gluon distribution function has no /i dependence, the rate for the process gg ggg 
will depend on \x through the Debye screening mass m 2 D which is used to avoid the divergence 
in both the scattering cross sections and the radiation amplitude and can be given by flTf 

Q ^2 

m l = ^T dpp[2f B + N f (f q + f 9 )] (10) 
n Jo 

In the above, Nf =2.5 is the dynamical quark flavour and g 2 = 4na s is the strong coupling 
constants. The above rate equations have been solved by using the distribution function as 
discussed below. 



B. Distribution function 

In order to evaluate energy, number densities, thermal quark mass and Debye screening 
mass we consider the approximations for the quark and anti-quark distribution functions 
given by 

f _ — \(g) e±T ^ \{o) e±X _ \ „±xreq / n \ 

~ \(5) e±/VT + eP/T ~ 1 + eP/T ~ m 
and gluon distribution function as 

u = \n q (12) 

where fj^L = (1 + e p ^ T )~ 1 and f e g q = (e P//T — For subsequent reference, we will call it 

as Modified Fermi Dirac type (MFD) approximation. It should be mentioned here that the 
most commonly used approximation (extended to non-zero p) is given by 
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The above approximation becomes Fermi Dirac distribution in the limit X q — > 1 and we 
will refer it as FD type approximation. It has been shown in appendix A that this FD 
approximation is not strictly applicable in case of an unsaturated plasma (i.e. for Aj < 1), 
though it has been used in all the earlier works [5]-|8|,|T3"|| for baryon free plasma (i.e. /i = 0). 



The MFD approximation (11) becomes FD type when \i — > 0. We would also like to compare 
the MFD approximation with the Boltzmann (BM) approximation given by 

/,(») = A 9 ( 9 -)e- ( ^ )/T (14) 

The BM approximation is found to be closer to the anti-quark distribution function at 
finite baryon density. Comparison between different approximations has been discussed in 
appendix A. 

We have evaluated the thermal quark mass m 2 q numerically using (9) for various approx- 
imations as well as for the Juttner distribution. Fig. 2 shows m? q /T 2 as a function of X q 
at two different values of x. Here, the gluon contribution to thermal mass has not been 
included since we are interested only in examining the difference arising from the use of 
different approximation for the quark distribution function. It may be seen that the FD 
approximation coincides with the Juttner distribution only in the limit X q — > 1. However, 
in the region of interest i.e. X q up to 0.1 or 0.2, the deviation is quite significant. Similarly, 
the Boltzmann approximation agrees with the Juttner distribution only at very small values 
of X q while MFD approximation agrees with the Juttner distribution up to X q « 0.35 even 
at x = 1.5. Since the plasma is highly unsaturated to begin with and the quark fugacity 
remains much below unity, we will use the MFD approximation in the subsequent calcu- 
lations. The advantage of using this approximation is that we can retain the factorisation 
used in [|] for the RHS of the rate equations (3) and (4) with the replacement of X q ( q ) with 
Aq(q) (see appendix B for detail). The rate equations are then solved in the following way. 

C. Formalism 

We assume the system to undergo a purely boost invariant expansion. Using the MFD 
approximation, (3) and (4) can be written as 

^("/) = ^ + - = VMl - \) - 2n 9 R 2 (l - ^) (15) 
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M„,un = + i = v?,(i - (i6) 

where R% = \o~ 2 n g and R3 = \cr^n g are the density weighted cross sections for the process 
99 ^ QQ an d 99 fi'fi'fi' respectively. Similarly, the equation for the conservation of energy 
and momentum can be written as 

1^ + ^ = (17) 
Or r 

where the viscosity effect has been neglected |2(J . In case of an ideal fluid, 3p —e. 



Using MFD approximation for quark, anti-quark (11) and gluon distribution functions 
(12), the energy density can be written as: 

e = T 4 [a 2 X g + b 2 (X q e x + X g e~ x )} (18) 

with a 2 = 87r 2 /15; b 2 = Nf(77c 2 /40) ; where Nf— 2.5 is the dynamical quark flavours. 
Similarly, the number densities for gluon, quark and anti-quark are: 



16 

n g = X g a l T z ] a x = — C(3) (19) 

7T Z 



n, = A,6iTV; h = ^ C(3)iV) (20) 

rig = A 9 -6iT 3 e- x (21) 

From the conservation of baryon number, one gets d^ijiBU^) = which results in 

ubt = (n q — riq)r = const (22) 

Using the expressions for the number and energy densities, the eq. (15,16) and (17,22) 
can now be written as 

^ + f + ! = «,(!- A.) -2H,(1 -^2) (23) 

Xg T T Xj 



^ + f +I = «,£^(l-^2) (24) 



A = A Q + ^ 7 ^ + ^(A Q -A ) (25) 
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3T 1 

1 T 



4 A 



\a 2 X g + b 2 X Q + b 2 Xq} 



(26) 



where 



A t = a 2 X g + b 2 (X Q + \ t 



The above equations are solved numerically estimating R 2 and R3 in a similar way as 
that of ref. |5|,[jj]. Following ref. f|, the rate R 2 is written as 



1 fiS 

i? 2 « .2AN f a 2 X g T ln(—) 

Ota A 



(27) 



where a s (= 0.3) is the strong coupling constant. The factor A = (A 9 + coshxAq/2) arises 
due to the thermal quark mass m 2 given by 



4na s T 2 
9 



[X g + X q cosh x/2] 



(28) 



which has been evaluated using MFD distribution functions. Following ref. , we estimate 
R 3 numerically 



R3/T 



32 a. 



N f 

X„ + Xg — - cosh a; 
6 



1 2 



3(2l Xg 

where I(X g , X q , x) is a function of X g , X q and x 



1 + 



2 m 



D 



9 T 2 



l(Xg,Xg,X) 



(29) 



J(A ff , A 9 ,x) 



dx 



s/Am 2 D 



dz- 



:i + z) 



cosh l (-Jx) 



cosh ^v^) 



[x + (1 + z)xo] 2 — 4:xzxn s - >K f J[l + x(l + -z)?/d] 2 — 4:Xzy D 



with ££) = m|) A/ and yo = mjj/s, s = 18 T 2 is the square of the average center of 
energy. The mean free path for elastic scattering is given by A/ 
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Aj 1 = - a s a x T X g 



N f 

Xg + Xg — - coshx 
6 



1 + 



2 m 



D 



9 T 2 



(30) 



mass 



(31) 



The above equations (29-31) are similar as that of used in ref. |7j except it has been rederived 
to have x dependence through the Debye screening mass m 2 D given by 

N f 



m 2 D = AnasT 2 



Xg + Xg — coshx 
6 



(32) 



derived using MFD distribution functions. It may be noted here that even in case of a baryon 
free plasma (x=0), the expression for m 2 D still differs from 47ia s T 2 X g which has been used 
in ref. The rates R 2 and R3 were calculated using the above equations. Fig. 3 shows 
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some typical results of the calculations of R2/T and R3/T as a function of X g for x=0.0, 1.0 
and 1.5. The values of X q = Xq were fixed at X g /5 . As can be seen, the equilibration rate R3 
initially increases and later on decreases with increasing baryon density whereas the rate R2 
is not affected much except at higher fugacities. This can be understood from the eq. (27) 
where R2 has a logarithmic x dependence whereas R% depends on x more directly through 
the Debye screening mass m 2 D (eq. 32). Also shown in the same figure are the results for 
R3/T used in ref. |F] for x=0 (dot-dashed curve) with m 2 D = 47ia s T 2 X g which does not 
include the quark and anti-quark contributions. It is seen that, with the inclusion of quark 
and anti-quark contribution in the Debye screening mass, the rates R3/T are lowered for 
X g > 0.2, but are enhanced for smaller X g values . 

D. Results 

The time dependence of X g , Xq, Xq and T were obtained by solving the set of rate eqs 
(23-26) numerically by fourth order Runga-Kutta method. We take the initial conditions 
from the Hijing calculations, T = 0.57 GeV, A 5 o=0.09, A g o =0.02 at r =0.31 fm as used in 
and treat xq as a parameter. We have carried out a parametric study to see the effect of 
chemical potential x on the chemical equilibration. 

Figs. 4(a-d) show the temperature T, chemical potential x — fi/T and the fugacities 
(X g , X q = Xq) as a function of r at few typical values of xq. Whether the plasma is baryon 
free or baryon rich, a common feature is that the process of chemical equilibration needs 
additional amount of energy which makes the plasma cool more rapidly than predicted by 
Bjorken's scaling T 3 r=const (dotted line, Fig. 4a). However, the plasma cools less rapidly 
in the presence of chemical potential compared to the baryon free plasma. Further, it is 
seen that the fugacity factors for gluon, quark and anti-quark do not reach the equilibrium 
values by the time the temperature T drops to the critical value T c « 0.2 GeV (Fig. 4(c) and 
(d)). It is important to note that the gluon equilibration rate slows down in the presence 
of finite baryon density. On the other hand, the quark and anti-quark equilibration rates 
increase slightly in the presence of finite baryon density. As shown in Fig. 4(b), the chemical 
potential x = fi/T also decreases with r. The overall effect, however, is that the plasma cools 
less rapidly in the presence of chemical potential as compared to the baryon free plasma. 

The different behaviour of X g and X q observed in Figs 4(c) and 4(d) as a function of r may 
be understood in the following way. The x dependence in the chemical equilibration arises 
due to the factors R 2 , R3 and baryon density (%) present in the rate equations (eq. 15 and 
eq. 16). As shown in Fig. 3, the rates R 2 and R 3 do not depend on x significantly. It was 



9 



found that even if R 2 and R3 are made independent of x, i.e. if we drop the term containing 
X q cosh x in m 2 , and m 2 D completely, the rate of chemical equilibration is not affected much. 
Therefore, the slowing down of the gluon equilibration rate is not due to R 2 and R 3 , but due 
to the dynamical evolution of the plasma which is affected in the presence of baryon density 
(See Appendix C for detailed discussion). The finite baryon density makes the cooling rate 
slower as the plasma evolves with an additional constraint of baryon number conservation. 
This also effects the plasma density which is mostly gluonic. It may be mentioned here 
that even though R 2 and -R3 do not depend strongly on x, the rate of chemical equilibration 
ultimately depends on how fast the net parton production approaches zero, i.e. how fast 
the RHS of eq. (15) or eq. (16) becomes zero. Since initially X g > \ q and -R3 > R 2 , a slight 
increase in gluon density produces more ggg as compared to qq or gg pairs ( i.e. more gain 
than loss in the gluon numbers as seen from eq. 15). As a result, the gluon equilibration rate 
decreases where as the quark or anti-quark equilibration rate practically remains unchanged 
until at a later time when it goes up slightly. The decrease in gluon equilibration rate will 
consume less energy and will make the plasma cool still slower. Therefore, the over all 
slowing down of the plasma cooling rate is due to both slower gluon equilibration as well as 
a slower expansion rate of the plasma that results in presence of baryon density. 

We have also investigated the rate of cooling by increasing the baryon density further, 
although our approximation breaks down at high x values. We find that the cooling rate 
never exceeds the Bjorken's limit which is consistent with the fact that for an ideal plasma 
Bjorken's scaling is the upper limit even in case of a baryon rich plasma. This may be 
verified directly from the conditions st= const and n j er=const (s and ris are entropy and 
baryon density respectively) without solving any rate equations. For an equilibrated plasma, 
x remains independent of r and the temperature follows the rule T 3 r = const. 

It is known that the rate of parton equilibration is enhanced considerably due to vari- 
ous factors like higher order gluon multiplication |pll , p2| , temperature dependent coupling 



constants |23| and viscosity. All these processes need to be included for a complete under- 
standing of the chemical equilibration during the hydrodynamical evolution of the plasma. 
We do not include all the above effects due to various complexities involved in the calcula- 
tions. For example, the perturbative calculations may not be valid if a s is allowed to vary 
with r which may have large value at the end of the evolution p3| |. Similarly, to calculate 
the rate for gg — > (m — 2)g; m > 5, one needs to understand the complete space time struc- 
ture of the multi gluon processes [^] . Similarly, inclusion of dissipative effects complicates 
the problem further. Therefore, in order to isolate the effect of baryo-chemical potential 
on the chemical equilibration , we have considered in this work a non-viscous plasma which 
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expands isentropically. We do not include any higher order gluon multiplication processes 
and also we use a constant value for a s =0.3. The results on the dilepton production yields 
are discussed in the following section. 



III. DILEPTON PRODUCTION 

Thermal photons and dileptons are the ideal probes to test the hydrodynamical evolution 
of a plasma created in heavy ion collisions at ultra- relativist ic energies. It was shown in ref. 
|14|JI7H that the thermal photon and dilepton rates are suppressed at finite baryon density 
when the rates are calculated at a fixed energy density. However, this suppression is primarily 
due to lower initial temperature that results at higher baryon density. We will show in 
the following that the presence of finite baryon density enhances the space-time integrated 
thermal yield of dileptons, in spite of the fact that the rate of production is suppressed |TJ 



. We have considered here the dilepton production, however this study can be extended to 
the thermal photons as well. 

The dileptons are produced predominantly via the reaction q + q~ — > . We ignore the 
annihilation and the Compton like reactions in the present calculations since their contri- 
butions may not be important for invariant masses above 1 GeV. The dilepton production 
rate dN/(d 4 xd A p) = dRjd A p (i.e. the number of dileptons produced per space time volume 
and four dimensional momentum volume) is given by: 

fp- = !Wf! ^fMfM** *'^ ~Vr- V2) (33) 

where v q q is the relative velocity between quark and anti-quark and a 1 ^ 1 is the total cross 
section for the reaction qq — > 



< = 3 W (34) 
M = p^p^ is the invariant mass of the dileptons. The above integral can be simplified to 



m 



dR 5 



/■OO />oo 

<IAP 24.i MV < M2 >i <MML m MM < 35 ' 

We evaluate the above integral numerically using the exact definition of quark and anti- 
quark distribution functions as given by (1). Fig. 5(a) shows the plot of dR/dM 2 as a 
function of dilepton invariant mass at fixed initial temperature T=0.57 GeV corresponding 
to an equilibrated (X q = 1.0) and a non-equilibrated (A g =0.2 and 0.02) plasma. As can 
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be seen, the rate of dilepton production is not sensitive to x when the plasma is highly 
unsaturated (the lower curve for \ q = 0.02). However, in case of an equilibrated plasma, 
it starts deviating more and more for higher x values (upper curve for A 9 =1.0). Since the 
plasma is highly unsaturated at the beginning (\ q m .02) and reaches up to \ q .2 as 
it equilibrates, the rate of dilepton production will be insensitive to the small values of 
baryo-chemical potential. It can be mentioned here that we have also estimated the thermal 
dilepton production rate using the approximation (11). We do not find much difference with 
the actual calculations particularly when the plasma is unsaturated. This further justifies 
the use of the MFD approximation both for quark and anti-quark distributions. Next we 
calculate the integral over the space and time using the expression: 



The rate dR/dM 2 is estimated numerically using the exact distribution functions. However, 
for the space time evolution of the plasma at finite baryon density (i.e. for T(r) and X q (r)), 
we use the results of the previous section which is obtained using the MFD approximation. 
Fig. 5(b) shows the above integrated yields at a fixed initial temperature of T = 0.57 GeV 
and initial fugacities \ g o = 0.09, \ q0 = \ q0 = 0.02 for different values of x . Although, the 
rate of dilepton production do not depend on x strongly, the space time integrated yield 
increases with x due to slower cooling of the plasma in the presence of baryon density . 
Since the energy density is directly related to the experimental observables, we have also 
calculated the dilepton yield at a fixed initial energy density of e = 9 GeV/fm 3 . This 
energy density corresponds to an initial temperature of T = 0.57 GeV at A 9 o=0.09 and 
Ago = 0.02 when x = 0. The results are shown in Figs. 6(a) and 6(b) for both the 
cases of equilibrated and non-equilibrated plasma. As shown in the figures, if the plasma is 
chemically non-equilibrated, the suppression of the dilepton yield at higher baryo-chemical 
potential is not significant even at higher invariant masses in contrast to the case for a 
chemically equilibrated plasma. 



In the present work, we have studied the effect of finite baryon density on the chemical 
equilibration of a longitudinally expanding quark gluon plasma. It is found that the rate 
of chemical equilibration for gluon slows down in the presence of finite baryon density in 
comparison to a baryon-free plasma, irrespective of various assumptions about the equili- 
bration mechanism . This results in a slower cooling of the plasma , which has an important 




(36) 



IV. CONCLUSION 
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consequence resulting in higher dilepton yields even though the rate of production may be 
suppressed at finite baryon density due to lower initial temperature. We have studied the 
thermal dilepton yields from the chemically non-equilibrated expanding plasma with finite 
baryon density . For calculating the dilepton yields, we have considered only an ideal non- 
viscous fluid that undergoes isentropic expansion and we have neglected higher order gluon 
processes as well as temperature dependence of the coupling constant. It is found that 
the space-time integrated yields of dileptons are enhanced if the initial temperature is held 
fixed. More importantly, even for a fixed initial energy density, the suppression of the yield 
at higher baryo-chemical potential is compensated to a large extent by the slow cooling of 
the plasma. 
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APPENDIX A 

Figs. A.l and A. 2 show the quark distribution function for p=0 and p=0.5 GeV as 
a function of X q at x—1.0. Figs. A. 3 and A. 4 show the corresponding plots for the anti- 
quark distributions. In all the cases, the FD approximations deviate significantly from the 
true Juttner distribution except when A — > 1. This deviation is more significant in case of 
the baryon free plasma although the FD approximation has been used widely in all earlier 
calculations. At x—0, the BM distribution (not shown in the figure), is found to be closer 
to the Juttner distribution at small fugacities, but it starts deviating even at small A when 
baryon density increases. On the other hand, the MFD approximation is relatively closer 
to the actual distribution over a wide range of fugacities (A.l and A. 2). However, the 
anti-quark distribution function (A. 3 and A. 4) becomes closer to BM approximation as x 
increases. At any momentum, for an unsaturated plasma at finite baryon density, the quark 
distribution function can be best described by the MFD approximation whereas the BM 
approximation is more suitable for the anti-quark distribution function. However, we do not 
want to use the MFD for f q and BM for f q as they will result in baryon asymmetry when 
x — > 0. Moreover, at finite baryon density, it is the f q distribution which dominates as the 
anti-quarks are strongly suppressed (A. 4). Therefore, we use the MFD approximation both 
for f q and f q in order to evaluate m 2 q1 m 2 D and quark and anti-quark energy and number 
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densities. We have also solved the rate equations using BM approximation both for quark 
and anti-quark distributions, but the final results do not change much . This is primarily 
due to the fact that the x dependence comes through the exponential e ±x in the numerator 
of eq. (1) and it does not matter if one uses BM or MFD approximations. However, we 
have used the MFD approximation for the calculations so that the same parametrization of 
ref. || can be employed by solving the rate equations. 

APPENDIX B 

For a baryon free plasma (/i = 0), Biro et. al. have used a simple factorisation for the 
rate equations. We can retain the same factorisation for both the processes gg qq and 
99 ^ 999 ( ec L-(3) and (4)) under the MFD approximations. Using the identity (7,8) and the 
MFD approximations (11,12) for quark, anti-quark and gluons, the gain and loss terms for 
the process gg # ggg (eq. (5) and (6)) can be combined to give 

r> _ P - fl 2 _ 1 \-U f rf 3 Pi f d 3 p2 r d s p 3 r d 3 p 4 v 

IX 9->9 1X Q^9 V A g Q Q) 2 J (2n)32E 1 ■> (2n) 3 2E 2 J (2n)32E 3 J (27r)32£ 4 

(2tt) 4 5 4 ( Pi +p 2 -p 3 -p 4 )E \M\ 2 ff(p 3 )fl q (p 4 )f^( Pl )f^(p2)e^ (B.l) 

where p\ + P2 = E\ + E 2 for mass less quarks or gluons. From the above expression, it is 
clear that the right hand side vanishes when X g = X q = Xq=l. Since the plasma is highly 
unsaturated, it may be reasonable to replace f eq with the Boltzmann distribution. The 
integral in eq. (B.l) can be written as 

I = \Swf J^^y^^fMfM (b.2) 

which represents the free space cross section for the process gg — > qq folded with the distri- 
butions for the initial particles. The cross section a gg ^ q q is given by 

It can be mentioned here that the above integral (B.2) is identical to the one which could 
have been obtained with the classical approximations, i.e. using the Boltzmann distribu- 
tion function for quark, anti-quark and gluon and eliminating the Pauli blocking and Bose 
enhancement factors in the final states in (5) and (6). Although, we use the quantum 
statistics, the same expression is obtained as the classical one due to the identity and the 
approximations eqs. (7,8,11,12). However Boltzmann approximation is used finally for f eq . 
Following H, eq. (B.2) can be factorised as 
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R 9 ^ 9 - R q -, 9 = \a 2 nj(l - ^) (B.4) 
Similarly for the rate gg — > ggg, one can use 

i? 2 ^ 3 - R 3 ^2 = ^3 nj(l - Xg) (B.5) 

where 

o- 2 =< a(gg -> >, tr 3 =< cr(## -> > 
V. APPENDIX C 

The quark gluon plasma in (1+1) dimenison expands following the scaling law T 3 r = 
const. If, chemical equilibration is not complete, the plasma cools faster than t -1 / 3 as 
additional energy is spent in chemical equilibration. Therefore, in case of a chemically 
equilibrating plasma, the cooling rate depends on the rate of equilibration of quarks and 
gluons present in the plasma. Further, as shown in the text, the presence of bar yo- chemical 
potential makes the gluon equilbration slow and also the temperature of the plasma falls 
slowly in comparison to the baryon free case. However, as will be shown here, the slow 
cooling of the plasma is not entirely due to the decrease of the gluon equilibration rate, 
but also due to the presence of baryo-chemical potential which effects the hydrodynamical 
expansion of the plasma. This has been demonstrated in fig. C where the temperature 
T, the gluon (A 9 ) and quark (X q ) fugacities have been plotted as a function of r both for 
baryon free and baryon rich cases (x—0 and 2.0, see the solid and dashed curves). As 
discussed in the text, the gluon equilibration rate slows down and the plasma cools slowly 
when the baryon density is finite where as the quark equilibration rate is not effected much. 
Next, we consider only a baryon free plasma, but we try to reduce the gluon and the quark 
equilibration rate by reducing the rate R3 and R 2 by a factor f=0.6 so that the plasma now 
cools with the same rate (see the dotted line) as that of baryon rich case with £=2.0 (dashed 
curve obtained using the normal value of R 2 and -R3). However, the equilibration rates for 
gluon and quarks now become much slower (dotted lines) as compared to the dashed curves. 
This indicates that if the plasma needs to cool with a rate as that of dashed curve just due 
to chemical equilibration alone , the gluon and quarks need to equilibrate with a rate much 
slower than what is shown in fig. C corresponding to the case of a baryon rich plasma. In 
other words, a higher equilibration rate of gluon and quark as that of dashed curves will 
consume more energy and the plasma will cool faster than what is shown by dashed curve 
(but still slower than the solid line) had the chemical equilibration been the only reason for 
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the deviation of the cooling of the plasma from the Bjorken's scaling. Therefore, in presence 
of baryon density, the slow cooling rate of the plasma arises due to two factors; one being 
the slow gluon equilibration rate which consumes less energy. The other one is due to the 
slow hydro-dynamical expansion as the plama now needs to expand conserving the baryon 
number as well. 
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Fig. 1 

(a) The product of quark and anti-quark distribution function as a function of p for 
different values of x at T = 0.57 GeV and 0.2 GeV. (b) Same as (a), but as a function 
of \ q . 

Fig. 2 

The thermal quark mass m 2 q /T 2 versus X q at x = 1.0 and x = 1.5 (the upper curves 
are increased by one unit). The solid curve is for the Juttner distribution function (eq. 
1 ). The other curves are obtained using different approximations, MFD (eq. 11), FD 
(eq. 13) and BM (eq. 14) respectively. 

Fig. 3 

The gluon production rate, R3/T, and quark production rate R2/T as function of X g 
(Aq = Xq = Xg/5) for x = 0.0, 1.0 and 1.5 . The dashed dot curve is for R3/T at 
x = 0.0 with m 2 D = 4iia s T 2 X g as used in ref. 

Fig. 4 

(a) The temperature, (b) the quark chemical potential x = fi/T, (c) the gluon fugacity 
X g , (d) the quark fugacity X q , as a function of r for x Q = 0.0, 1.0 and 1.5 with the 
initial conditions T = 0.57 GeV, X gQ = 0.09 and X q0 = 0.02. The dotted line in (a) 
corresponds to the temperature as per the Bjorken's scaling. 

Fig. 5 

(a) The dilepton production rate dR/dM 2 as a function of invariant mass M for X q = 
0.02, 0.2 and 1.0 . The initial temperature of the plasma is To = 0.57 GeV. (b) The 
space time integrated dilepton yield versus invariant mass M for a chemically non- 
equilibrated plasma (A 90 =0.09, A g0 =0.02 ) at a fixed initial temperature of T = 0.57 
GeV with x = 0.0, 1.0 and 1.5. 

Fig. 6 

(a) The integrated dilepton yield from a chemically equilibrated plasma at a fixed 
initial energy density for various values of xq. (b) Same as Fig. 6(a), but from a 
chemically non-equilibrated plasma. 

Fig. A 

The quark (f q ) and anti-quark (f q ) distribution functions versus X q at p = (A.l and 
A.3) and p = 0.5 GeV (A.2 and A.4). 
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